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The general expression of the angular distance between two point sources as measured by an 
arbitrary observer is given. The modelling presented here is rigorous, covariant and valid in any 
space-time. The sources of light may be located at a finite distance from the observer. The aberration 
and the gravitational deflection of light are treated in a unified way. Assuming the gravitational 
field to be weak, an explicit expansion of the angular separation within the post-post-Minkowskian 
approximation is carried out. The angular separation within the post-Newtonian approximation 
truncated at the order 1/c 3 is straightforwardly derived. 
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I. INTRODUCTION 



An accuracy of the order of the microarcsecond (/xas) is expected in future space astrometry missions, like the Global 
Astrometric Interferometer for Astrophysics (GAIA) which is planned to be launched at the latest in 2012 ffl, IH. fiol] . 
f**«> \ or the Space Interferometry Mission (SIM) [111] . Several models of positional observations compatible with such a 
level of accuracy have been proposed within the post-Newtonian approximation of general relativity (see, e.g., [(| and 
[H). An even more general model was formulated by Klioner 0] in the framework of the parametrized post-Newtonian 
(PPN) formalism [14j |. with parameters (3 and 7. Such models are sufficient for the currently envisaged astrometry 
space missions. However, quite ambitious laser interferometric missions like LATOR [13] and ASTROD are now 
proposed with the aim of testing the post-linear approximation of metric theories near the Sun. The main purpose of 
the present work is to outline a rigorous and general modelling of the angular distance between two light sources as 
1 , measured by an observer having a given motion and then to specialize to the different approximations which will be 
H ■ needed in the foreseeable future. 

All the perturbation expansions performed in this paper are deduced from the rigorously covariant formulae given in 
Sec. II. It must be emphasized that the aberration and the gravitational deflection of light are not separately treated. 
Moreover, the sources of light (supposed to be pointlike) may be located at finite distances from the observer. As a 
consequence, all the formulae given here take into account the motion of the observer and are suitable both for remote 
C3 ' stars or quasars and for objects moving in the solar system. 

In Sec. II we derive the general expression of the angular separation between two point light sources as measured by 
a given observer. We use in this section a coordinate-free approach totally avoiding the index notation owing to the 
simplicity of the calculations. In Sec. Ill we show that the angular separation can be explicitly calculated when the 
reception time transfer function as defined in [8( is determined provided that the motion of the observer is known. In 
Sec. IV, we get the general expansion allowing to determine the angular separation within the post-post-Minkowskian 
approximation. In Sec. V, we carry out the calculations within the post-Newtonian approximation. Section VI 
delivers some concluding remarks. 

Throughout this work, c is the speed of light in a vacuum and G is the Newtonian gravitational constant. The 

Lorentzian metric of space-time P4 is denoted by g. We adopt the signature (H ). As long as explicit coordinates 

are not needed, four dimensional vectors are denoted by usual letters : u,v,l,l' ,U, f3 u / v , f3 v / u , .... We put u.v = 
g(u,v) — g a pu a vP and u 2 = u.u. 

When an explicit introduction of local coordinates is required, we suppose that space-time is covered by some 
global coordinate system x a — (x°, x), with x° — ct and x — (x l ). Moreover, we assume that the curves of equations 
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x % = const are timelike in the neighbourhood of the observer. This condition means that goo > in the vicinity 
of the observer. We employ the usual three-dimensional vector notation a in order to denote either the ordered set 
(a 1 , a 2 , a 3 ), or the orderer set (01, 0.2, (Z3). Since there cannot exist any confusion with the above-mentioned notation 
of the scalar product of two 4- vectors, given a = (a 1 , a 2 , a 3 ), for instance, we denote by a.b the quantity a l b l if 
b = (6 1 , 6 2 , 6 3 ) and a l bi if b = (61,62,63)), the Einstein convention of summation on repeated indices being used 
in each case. The quantity \a\ denotes the ordinary Euclidean norm of a : \a\ = (SijO 1 ^) 1 ' 2 if a — (a 1 , a 2 , a 3 ), 
and \a\ — (<5 y a^et.,) 1 / 2 if a — (01,02,03). The indices between parentheses characterize the order of a term in a 
perturbative expansion. Theses indices are set up or down, depending on the convenience. 

When there exists a risk of ambiguity, we use the subscript o for quantities related to the observer : for example, 
x Q is the point-event where the observation is located. We use the subscript e for quantities related to the emission 
of a light ray : t e is the coordinate time of emission and x e denotes the position of the emitter. 



II. ANGULAR DISTANCE AS MEASURED BY A GIVEN OBSERVER 

Let us consider a light ray T received at point x a = (ct ,x ) and let us recall how is defined the direction of this 
ray as measured by an observer 0(u) moving at x Q with a unit 4-velocity u. The three-space relative to the observer 

0{u) at point x is the subspace H Xo (u) of tangent vectors orthogonal to u. An arbitrary vector V at x admits one 
and only one decomposition of the form 

V = V\ U + V ±u , (1) 

where Vji is colinear to the unit vector u and V± u is a vector of the three-space H Xo (u). Since V± u and u are 
orthogonal, one has 

Vj|„ = (u.V)u (2) 

and 

V ±u = V - (u.V)u . (3) 

The vector V± u is called the (orthogonal) projection of V onto the three-space relative to the observer 0(u). Its 
magnitude \V± U \= y / —V± ll .V± u is given by 



|Vi„h ^{u.V) 2 -V 2 . (4) 
The direction of vector V as seen by the observer 0(u) is the direction of the unit spacelike vector V± defined as 

V * Y± - v -( u - v ) u 

x « \V±.\ ^{u.V) 2 - V 2 ' {> 

Consider now a light ray T received at x a and denote by I a vector tangent to T at x . In this work, we always 
assume that a vector tangent to a light ray is a null, future directed vector, so that 

Z 2 = 0, u.l>0. (6) 

The direction of the ray T as measured by the observer 0(u) is the direction of the vector V]_ . By using Eq. ([5]), 
and then taking into account Eqs. ((6|), it is easy to see that 

'i« = i— (7) 

Let r' be another light ray received at x Q . If I' denotes a vector tangent to T' at x , the direction of V as observed 
by 0(u) is given by Eq. (|7|) in which I' is substituted for I. As a consequence, the angular separation between T and 
r' as measured by 0(u) may be defined as the angle 4> u between the two vectors l\ and l± belonging to the same 

subspace ni 3 J (w) (see, e.g., @,[l3). The angle <p u may be characterized without ambiguity by relations as follow 

cos <f> u = -l* x .l'L, O<0 U <7T. (8) 
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By using Eq. {7) and then taking into account that it is a unit timelike vector, it is easily deduced from ([SJ that 

l.V 

cos (p u = l- ? ,w , A , (9) 
(u.l)(u.l') 

As a consequence, the angle <p u is determined by 

Si 44MM' 0<<t> u <n. (10) 



Since Z and V are supposed to be null vectors, we have 



U' = -\(l'-lf. (11) 
By substituting for l.V from Eq. (jTTJ) into Eq. (fT0|) . we get the formula 

sin 2 ^ - 1 (*' ~ 2 (12) 

2 " 4(u.l)(u.V)' ' 

Let us now compare <p u with the angular distance <p v between T and T' as measured by an observer moving with a 
unit 4-velocity v. It follows from (JXDJ) that 

2 (u.l)(u.V) 2 v ; 

The factor ahead of sin 2 4f is the aberration factor due to the motion of the observer 0(u) relative to the observer 
0{v ). Let us determine this factor in a more explicit form. 

Denote by r the proper time of the observer <D(u). The spatial displacement of 0(u) in the 3-space (v) relative 
to 0(v) between instants r and r + dr is the vector 

dx± v — c [u — (u.v)v\ dr . 

This displacement corresponds to a timelike displacement in a frame comoving with 0(v) given by 

dx\\ v = c{u.v) vdr . 

As a consequence, the velocity vector w u / v of the observer 0(u) relative to the observer 0(v) may be defined as 

Wu/v = cPu/v , (14) 

where (3 U / V is the vector given by (see also 0) 

dx ±v u 

^=1^1=— v~ V - (15) 
The magnitude of the spacelike vector (3 U / V has an expression as follows 



(16) 



This last relation is equivalent to 

-r^ = i- \Pu h \ 2 • (17) 

So, u.v is the 7 factor corresponding to the velocity Iwn/ul- 

Conversely, the velocity vector of the observer 0(v) relative to the observer 0(u) is given by 

v 

Wv/u = C(3 V / U , Pv/u = u . (18) 
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It results from Eq. (fTBj) that \/3 u / v ) \ and \/3 v /u I have the same magnitude : 

\(3 v /uHPu/v\, (19) 

a property which is well known in special relativity. Let us note, however, that (3 U / V and (3 V / U cannot be opposite 
vectors since they do not belong to the same subspace, except in the case where u = v at point x a . 

We are now in a position to write Eq. I|13p in a more explicit form involving [3 U / V and the unit spacelike vectors lj_ 
and which characterize the directions of the rays T and r", respectively, in the three-space relative to the observer 
0(v). Indeed, I and I' may be written as 

l = (v.l)(v + l* ±v ), I' = (v.l')(v + l'lj (20) 

and Eq. (fl~5f yields 

u= (u.v) (v + (3 u/v ). (21) 

Since lj_ , l'± and j3 u / v are orthogonal to v, Eqs. (|20p and (f2"Tj) imply 

u.l = (u.v)(v.l){l+p u/v .l* ± J, u.l'=(u.v)(v.l')(l + u/v .l'l v ). (22) 

Substituting for u.l and u.V from their respective expressions given by Eqs. (|22|) into Eq. (TT3|) and taking Eq. (fT7|) 
into account yields the fundamental relation 

. 2 0u 1_ I 2 • 2 0u / oq n 

sin — = t r-; - — r- sm — , (23) 

which gives (/)„ as a function of 0„ and of vectors which are all in the three-space Tl^] (v) relative to the observer O v . 

Since the null vectors I and V are assumed to be future directed, the directions in which the observer 0(v) is seeing 
the sources of the light rays V and V are characterized by the unit vectors N v and N' v defined as 

N v = - l* ±v , K = ~ l 'L ■ ( 24 ) 
Of course, <j> v is the angle between N v and N' v . 

Denote by ip v (resp. ip' v ) the angle between (3 U / V and N v (resp. f3 u / v and A^). One has 

\(3 u/v | cosip v = -N v .(3 u/v = f3 u / v .l* ±v , 0<ip v <tt, (25) 
| /3 U/V | cos = -K.j3 u/V = u/v .l'l v , < tf v < n . (26) 

Then Eq. (12"3")) may be rewritten as 

2 (1+l/^lcos^Xl+IAvJcosVO 2 ' 

Conversely, put 

tf« = K = -i'Z t , (28) 

and define the angles and by 

\Pv/u\ cosV>« = -N u .(3 v/u = ^ v / u .l± u , 0<ip u <n, (29) 

\P v/u | COS< = = Pv/u-l'L , < < < 7T. (30) 

We may write 

. 2 1_ IA>/« I 2 • 2 <K /oi\ 

sm — = ; ; r-r : : r Sm — . 31 

2 (1+ \I3 V/U \ cosVv)(l+ \{3 v /u\ cosiP'J 2 
It is worthy of note that comparing Eqs. (f2"T)) and (f3"Tj) yields an identity as follows : 

(1+ \/3 v /u \ cos-0 M )(l+ \/3 v/u \ cost/4) 1- |Ay„| 2 



1- \Pv/u\ 2 (1+ |At/«l COS<0„)(1+ \Pu/ v \ COS1p' v ) 



(32) 
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III. ANGULAR SEPARATION EXPRESSED IN A GIVEN COORDINATE SYSTEM 



Up to now, we have not specified any coordinate system. Henceforth, we shall assume that the coordinate system 
x a = (x°, x l ) is chosen so that d = d/dx° is a timelike vector in the neighbourhood of the point of observation. This 
condition means that g o > in the vicinity of x = (x®,x l ). On this assumption, we can consider that the curve 
defined by x l = x % = const is the worldline of an observer 0(U) at rest relative to the reference system defined by 
the coordinate system x a and passing through x Q . At x a , this observer is moving with a unit 4- velocity U given by 



U 



1 



a 



(33) 



V5oo dx° ^/goo 
A null vector cannot be orthogonal to a timelike vector. As a consequence, the null vectors I and I' are such that 

l = (l.d )^0, l' = (l'.d )^0. (34) 
So the above formulae remain valid under the substitutions 

I — ► t, I' 

where I and I' are the null vectors defined by 

I 



?. 



(35) 



I 



(l.do) 



I' 



I' 



(I'.do) 



(36) 



The introduction of vectors I and V will be justified below. 

Using these definitions, the angular separation (f>u between T and V as measured by the observer 0(U) is given by 
a relation as follows 



sin 



2 <PU 



1 {I'-lf 



and the relation between </> u and <pu may be written as 

. 2 k (uT)(uT') . 2 </>u 

sin — = — ^— sin 

2 {u.l){u.V) 2 

Noting that the covariant components of I and V are given by 



l' n 



V 

11' = -!- 



and then, taking into account Eq. 



it is easy to see that Eq. (|3T|) reads 
Then, applying the general formula ([?7)) with v = U yields 



. 2 VU 

sm 

2 



sm 



I A 



IV I 



4(1+ | Pu/U I COS Ipu) (1+ I u /u I COS'i/'y) 



w j «-Ji)(i;--ii 



(37) 



(38) 



(39) 



(40) 



(41) 



In many practical problems, however, it will be indispensable to express all the factors involved in sin ^ in terms 
of coordinate-dependent quantities. Of course, the contravariant components of the vector (3 u /u could be explicited 
by replacing v by U in Eq. (Til)]) , and then using Eq. (|3"3")l . Nevertheless, using Eqs. (f3"3"| and ([3"5j) - (T4TH) yields 
straightforwardly 



sm 



1 g^^-^)^.-; 3 ) 
K) 2 (i+^r TO )(i+^f;) 



(42) 
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where quantities f3 l are the components of the coordinate- velocity of 0(u) divided by c : 

dx l dx % 



dx° cdt 



Since u is a unit timelike vector, one has 



g 00 + 2g 0l (3 l + g %] fi l ^ 



Finally, substituting for l/(u ) 2 from Eq. (j4"4")l into Eq. (j4"2")l yields the fundamental formula 



sin 



2 Vu 



(goo + 2g ofc /? fc + g k iP k P l ) g^ti - U){j' 3 ~ lj) 



(43) 



(44) 



(45) 



This formula gives the expression <f> u as a function of the metric, the coordinate velocity of the observer 0(u) and 
the six quantities li and l[ at point x a . In order to determine the quantities li, we may proceed as follows. Suppose 
that the light ray T received at point (ci , x ) is emitted at point (ct e , x e ). The travel time t Q — t e of this ray may be 
considered as a function of x e , t Q and x a . So we can put 

t Q t e Ti*(x e: £ , *^o) : (^^0 

where T r may be called the reception time transfer function. Then, as it is proved in Q , we may write 



dx i 



1 



dto 



(47) 



As a consequence, once the reception time transfer function is explicitly known in a given space-time, it is possible 
to obtain the angular separation as a function of the coordinates of the emission and reception points, and of the 
coordinate speed of the observer. 

Henceforth, we focus our attention on the case of a weak gravitational field. So we suppose that space-time is 
covered by some global quasi-Galilean coordinate system x a — (x , x) in which the metric may be written as 



gp,v = r)i^ + , r)^ = diag (1, -1, -1, -1) . 



(48) 



where the gravitational perturbations meet the conditions \h^\ << 1. According to this assumption, the time 
transfer function T r may be written in the form 



1 



^r(*^e; Xq} — \x e X Q \ -f- Q r (x e: t Ql IC ) 



Substituting for T r from Eq. (|4*9")l into (|4T|) and then putting 



j\P = -2- N H 



it is easily seen that 

(li) = N i + X i (x e ,t ,x ), (ity = N' i + X i (x' e ,t ,x ) 
where the functions Xi(x, t , x ) are defined as 



1 - 



d<d r (x,t ,x ) 
dto 



dQ T {x : t Q: Xo) dQ r (x, t : *^o) X Q 



dxl 



dt 



(49) 



(50) 



(51) 



(52) 



The functions Xi(x, t a , x D ) characterize the fact that a light ray emitted at point x e is not arriving at (ct ,x in 
the direction defined by the direction cosines {x l — x a )/ \ x e — x a |. So we call them the (gravitational) deflection 
functions. 
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For the sake of brevity, we introduce the vector notations 

i={k}, i' = &}, (53) 



N = m = JHL_T2L N > = {N >i }= A_^ } (54 ) 

|«Ce *Eo\ l*^e *Eo\ 

A = A(a: e , i D , x Q ) = {\i(x e , t D , x )} , A' = \(x' e , t , x D ) — {\i(x' el t , x D )} . (55) 

Thus Eqs. (|5ip may be written in the form 

l = N + X, 1' = N' + A' . (56) 

It is clear that N is the zeroth-order, unit coordinate vector from the observer to the source. Noting that |JV'— 7V| = 
y/2(l — N.N'), we define a "unit vector" s characterizing the angular separation as 

N' - N 

s = {s 1 } = - . (57) 

y/2(l- N.N') 

Finally, we put 

/3 = {/?}, /3 2 =%/?^. (58) 

Substituting for li and l[ from Eqs. (|5ip into Eq. (|37p . and then using the above-mentioned notations, we get 

. 2 4>u 1 1 - P 2 + h 00 + 2h ok (3 k + h kl (3 k (3 l 
sin 



2 2 [1 + /3.(N + A)][l + [3.{N' + A')] 

x |(1 - N.N') (l - jfe ij sV') + ^2(1 - N.N') [(A' - A).s - fc^(A- - Xi)s j ] 



+-(A' - A) 2 - -fc^ - A0(A; - Xj) , (59) 



2 V ; 2 

where the quantities W are defined by 

g ^ = ^ + . (60) 
Let us emphasize that the formula (159|) is rigorous. 

IV. ANGULAR DISTANCE WITHIN THE POST-POST-MINKOWSKIAN APPROXIMATION 

Let us assume now that the metric perturbations may be represented at any point a; by a series in ascending 
powers of the Newtonian gravitational constant G 

oo 

h lw {x,G) = Y i G n g^){x). (61) 
i 

We shall put 

h$=G»g${x). (62) 
More generally, given any function f(x, G) which admits an expansion as follows 

oo 

f{x,G) = Y t G n vF>{x), (63) 

n=0 

we shall put for the sake of brevity 

f {n} {x,G) = G n u M {x), (64) 
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so that / may be written as 



f(x,G) = J2f 



in) 



(65) 



n=0 



For the sake of convenience, we shall sometimes use the notation fr n ) instead of /("). 

In what follows, we restrict our attention to the post-post-Minkowskian approximation, which means that we neglect 



all terms involving hffi , h^), ... As a consequence k lJ may be replaced by 



jfe« = kf 1} + fc« + 0(G 3 ) . 



where 



uij _ _ 1,(2) 1.(1)1,(1) 



(66) 

(67) 
(68) 



According to Eq. (|61[) . we assume that the vector deflection functions A and A' admit expansions as follows 

A = A (1) +A (2) +0(G 3 ), A' = A' (1) +A' (2) +0(G 3 ). (69) 

Furthermore, we point out that in all practical cases, the speed of the observer is much smaller than the celerity of 
light in a vacuum. So we henceforth assume that 



which implies that a condition as follows 



< /3 2 < 1 



(1 + /3JV)(1+/3JV')^0 



(70) 



(71) 



is satisfied. 

Using the previous expansions, a long but straightforward calculation shows that the angle cj> is determined by an 
expansion as follows 



sin 



l-(3 2 



1 



2 (1+/3JV)(1+/3JV') 



- N.N' ) 



V2(l- JV.JV') + £ (2) ] + i (A'« - A«) 2 + 0(G 3 ) 



where 



- 41 j = f3^2 + h V s 3 

B« = (a'M-aw).*, 



1+/3JV 1 + /3JV' 



(72) 

(73) 
(74) 



and 



^ 2 ) = 



J$ + 24 2) /3 l + fcg^ + 2h$p i + hVf?? ( p.\& /3.A'« 



1-/3 2 1~(3 2 \l+f3.N ' 1+p.N' 

'h$ + 2h$p k + h$0 k l /3.A« /3.A'« 



l~f3 2 



p.X^ /3.A'( 2 ) 



/3.A« 



1 + p.N 1 + 0.JV' J ij V *t J ' 5 
/3.A« /3.A'« / /3.A'« 



1 + p.N 1+p.N' \l+p.NJ 1+p.N 1+p.N' \l+p.N> 



(75) 



B (2) 



y(2) _ A (2) 



+ 2h$/3 k + h$f3 k f3 l p.\W p.X'W 



1-/3 2 



1 + P.N 1+p.N' 



A'(D _ A (D 



(i) _ 



(76) 
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Denote by 4>u^ the angle denned by 



6 (0) 

2 Yu 



1 



l-(i 2 



2 2 (1 + p.N)(l + p.N') 



-(1-N.N'), 0<^°)<tt. 



(77) 



It is clear that </4°^ is the angle which would be measured by the observer 0(u) in the absence of any gravitational 

(78) 



field. As a consequence it may be assumed that <f> admits an expansion as follows 

0„ = 0(0) + 0« + 0(f) + 0(G 3 ) . 
Using Taylor's theorem, it may be seen that 

d°) 



■ 2 Yu .2 

sm — = sm 
2 



±- + \ 0« sin^> + \ rf£0 sin0i°) + \ (^) 2 cos^ + 0(G 3 ) . (79) 



Now compare Eq. (|79| with Eq. (|42|) . and then take Eqs. ((73|) - ([76|) into account. We find for the first-order 
perturbation term 



2 v^I+^jvxT+^ivo 



(80) 



or, more explicitly 



<t>V = tan%- 



h® + 2h$/3 i + h^;r:p 



(1) ^ p.\W /3.A'« 



l-f3 2 



h-j's s> - 



1 + p.N 1+p.N' 



yj{l + P.N)(l + P.N') 

For the second-order term, we obtain 



cos ■ 



A'W-AW ).s. 



(81) 



AO) 



2 vPPOfiW) 



p(2) 



1 



1-/3 2 



2(1+/3JV)(1+/3JV') sin ^ 



(0) 



,(0)' 



A'W-AW) - h-tan 2 ^-] f A'« - A« ) . s 



(82) 



where 



do)' 



= A ™ - ~ 1 - tan 2 ^- ) U W 



,(o)' 



1 -/? 2 



- 1 - tan z 



1-/3 2 



-- [ 1 + tan 1 



\ h$ + 2/$ ) /3 i + / /3.A« /3.A'« 



1-/3 2 



1 + P.N 1 + p.N' 



,(2) 



1 



1 + tan z 



« 2 

/3-A( 2 ) P.X'^ , I , . 

— 3 + tan 

1+P.N 1 + p.N 4 \ 2 



<^>\ ( h^ + 2h^ + h^ll_ p.xw p.yw \ (1) _ (1) (1) 

2 M 1 - /3 2 1 + p.N 1 + p.N j v V ^ ^ 



6 (oy 



/3. 

1 + p.NJ + \l+p.N> 



1 



2 <ru 



-ll + tan —l 1+/3JV - 1+/3JV ,, 



(83) 
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and 



p(2) 



B 



(2) 



1 - tan' 



(A'( 2 )-A( 2 )). S + i (\ + tan 2 ^ 



'h$ + 2h$0 k + h$0 k l /3.A« /3.A'« 



'l-tan 2 ^) hfsW 



l-(3 2 
(a'« - \V 



1 + p.N 1 + /3JV' 



(84) 



Post-Minkowskian approximation. The previous formulae drastically simplify if we neglect all the contributions of 
second order with respect to G. Let us introduce the unperturbed angular distance <py which would be measured by 
the observer 0(U) at rest relative to the chosen coordinate system in the absence of gravitational field : 



6 {0) 

2 



1 



Mir — = -(1 -N.N'), O<0^<tt. 



We obtain within the linearized, weak field approximation 
1 1 - P 2 



sin 



2 Yu 



2 2 {1+/3.N)(1+/3.N' 

p.x p.y 



-1(1 -N.N') 



h 00 + 2h 0l [3* + hi,?? . i j 

1 + T^J 2 ijS 



1 + P.N 1 + p.N' 



v/2(l - N.N') (A' - A), s + 0(G 2 ) J . 



Assuming (fry ^ and \/3\ << 1 allow to develop the square root of the r.h.s. of ([86]) as follows 
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21 + p.N 21 + p.N' 



l(X'~X).s + 0(G 2 ] 



(87) 



V. ANGULAR DISTANCE WITHIN THE POST-NEWTONIAN APPROXIMATION 



Let us now suppose that the weak gravitational field considered in the previous section is generated by a system 
of N sclf-gravitating bodies within the slow-motion, post-Newtonian approximation. By choosing a convenient quasi 
Cartesian coordinate system (e.g., a standard post-Newtonian gauge), the metric may be written in the form 

2w 

hoo = T + 0(c- 4 ) , 

-3\ 



hot = 0(c~ 6 ) 



w 



2 7 -^+0(c- 4 ), 



(89) 
(90) 



where w is given by 



w(x°,x) = G 



N 



pa ^°^d 3 x' + 0(c- 2 ) 



\x — X' 



(91) 



p a being the rest mass density of the body a and T> a the domain occupied by the body a at instant x°. 

Within this approximation, we may neglect all the terms which are of the order of 1/c 4 and we have to retain only 
the terms of order 1/c 2 in the gravitational deflection functions A^ and A^. Substituting for h^v from Eqs. 
into ([86]) yields 
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(92) 



11 



Hence, assuming (j>ffl ^ and taking the square root of the r.h.s. of Eq. ([92 



d°) r 



2 ^/(l + (3.N)(l + (3.N>) 



l-(l + 7) 



w 1 



/3.(A + A') 



(A' - A) . s + 0{c 



Using the results established in Sec. IV, it may be seen that the angle <f> u admits the expansion 

K = 0W> + + 0(4) , 

where <^ is defined by (|77|) and (jJh N ^ is given by 

/3.(N + N') 



1 



1 + N.N' 



2 ,w ,v „,„ ,w / 1 — N.N' 
7 (A' -A). a -2(1 + 7)- 



1 + AT. AT' 



1 + N.N' 



1 - AT. AT 



1 + AT. AT 
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(93) 



(94) 



(95) 



This formula yields a thorough determination of the influence of the motion of the observer and of the gravitational 
light deflection on the angular distance described within the post-Newtonian approximation. 



VI. CONCLUDING REMARKS 



This work yields a general, rigorously covariant determination of the angular distance between two point sources of 
light as measured by an observer moving with a given velocity. The fundamental formulae (|45[) and (|47|) show that the 
theoretical calculation of the angular distance can be carried out when the so-called reception time transfer function 
is known throughout space-time. So it can be said that the problem of space astrometry involves one and only one 
unknown function. Furthermore, it must be pointed out that our formulae do not involve any tetrad formalism to 
define the directions of light rays. It is also worthy of note that the aberration and the gravitational deflection of light 
are not separately tackled. 

The case of weak gravitational fields is explicitly and thoroughly treated within the post-linear approximation. The 
formulae (|78p - (l84)) yield the perturbative expansion of the angular separation up to the second order in G. These 
results will enable to discuss the effects which could be observed by the next generation of stringent tests of general 
relativity like LATOR or ASTROD. Finally, the expression of the angular separation is obtained within the post- 
Newtonian approximation truncated at the order 1/c 3 . We are currently using this last calculation for determining 
the contributions of the mass-multipole coefficients of the gravitational field to the bending of light in the vicinity of 
a giant planet like Jupiter. 
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